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Abstract 

Confined and deconfined phases are defined tlirough nonperturbative correlators and 
nonperturbative background perturbation tlieory is used to compute tlie critical temper- 
ature Tc and spatial string tension. Taking evolution along one of the spatial axes the set 
of Hamiltonians H{ni,n2) with different Matsubara frequencies is obtained. 

Meson and glueball screening masses and wave functions are computed for H{0,0) and 
compared with available lattice data. 



1. Finite temperature QCD, and in particular the dynamics of deconfinement phase transi- 
tion and the physics beyond the critical temperature is an important testing ground for our 
ideas about the nature of the QCD vacuum. 

During the last few years, a lot of lattice data have appeared which point to the nonpertur- 
bative character of dynamics above T^. To this belong i) the area law of spatial Wilson loops 
[1,2], ii)the screening masses of mesons and baryons [3-5] and glueballs [6], iii)the temperature 
^ - dependence of Polyakov-line correlators [7-9]. In addition, the behavior of e — 3p above 
.has a bump incompatible with the simple quark-gluon gas picture [10]. Thus the inclusion of 
'nonperturbative (NP) configurations into QCD at T > and also at T > Tc is necessary. 

One of the authors has developed a systematic method for QCD, treating NP fields as a 
background and performing perturbative expansion around that both for T = [11] and T > 
[12]. 

To describe the phase transition, a simple choice of deconfined phase was suggested [12] 
where all NP color magnetic configurations are kept intact as in the confined phase, whereas 
color electric correlators that are responsible for confinement vanish. 

This picture together with the background perturbation theory [11,12] forms a basis of 
quantitative calculations, where field correlators (condensates) are used as the NP input. 

As an immediate check one computes for a different number of flavours in good 
agreement with lattice data. One purpose of the present letter is to treat quantitatively the 
points i)-iii) systematically in the framework of the same method. Doing so we meet some 
limitations. First, in the description of the deconfined phase in the spirit of the aforementioned 
method [12], one uses the magnetic correlators ( and spatial string tension as) from the lattice 
calculations [1,2]. In the confined phase below Tc and up to the temperatures T < 1.5Tc in the 
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deconfined phase these correlators are practically the same as are known for T = 0. Therefore 
in the mentioned temperature region the spatial string tension Us is regarded as T- independent 
and coincides with the standard string tension a at T = . This is because the main dynamics 
(and contribution to the free energy) is due to the scale anomaly which essentially changes at 
the temperature of the order of the dilaton mass, i.e. at T < IGeV. However, at large T, a new 
physics sets in: that of dimensional reduction, which causes a.^ to grow [2,13]. For temperatures 
larger than 1.5 we have used the scaling behavior of o"s(T) from [2], that is consistent with 
the dimensional reduction regime. We confine ourselves to the region < T < 2 ST^ and 
consider problems i)-iii) point by point. 

The second limitation is of technical character, since the proposed method is used for sim- 
plicity in conjunction with the 1/A^c expansion. Therefore values of are computed in the 
leading and sub leading orders of 1/Nc expansion and strictly speaking the conclusion obtained 
for the first order of deconfinement phase transition cannot be extended to Nc = 2. 

To calculate the screening masses one should first write the Green function of quark- 
antiquark (qq) or two gluons {gg) in terms of the Wilson loop average [14,15] and derive from 
it the Hamiltonian taking evolution along one of the spatial axes. In fact, one obtains a set of 
Hamiltonians, containing all higher Matsubara frequencies. Our results for mesons in general 
agree with one of the choices adopted before in [17], and justify in that special case the confining 
ansatz made there. 

Of special importance are the gg Hamiltonian and the corresponding glueball screening 
masses (GSM). Here the lowest Matsubara frequency is zero and the next one is 27rT. Therefore 
within our picture only the nonperturbative interaction between gluons - spatial string tension 
- determine the temperature behavior of the lowest GSM. In contrast to that, the meson 
screening masses for the lowest states have an additional temperature dependence through the 
lowest Matsubara frequency. 

We calculate the GSM corresponding to lowest states and compare our results with recent 
lattice data [6]. With the scaling behavior crs(T) [2] GSM forms the similar pattern to that of 
the meson case. This picture is in agreement with existing lattice data [7-9]. 

2. We derive here the basic formulas for the partition function, free energy and Green 
function in the nonperturbative background formalism at T > [12]. The total gluonic field 

is split into a perturbative part and a nonperturbative background 5^ 

A^ = B^ + (1) 

where both 5^ and are subject to periodic boundary conditions. The principle of this 
separation is immaterial for our purposes here, and one can average over fields and 
independently using the 'tHooft's identity^ (for more discussion see [18]) 



/ 



/ DB^r]{B) JDa^exp{- S{B + a)) 
JDB.r^iB) 



DA,exp{-S{A)) = J--,'iy-JJ--,-J^--y-^'^JJ (2) 



=<< exp{—S{B + a) >a>B 

with arbitrary weight r]{B). In our case we choose r]{B) to fix the field correlators and the 
string tension at their observed values. For simplicity have omitted gauge fixing and ghost 
terms in (2). 

To the lowest order in gaf^ the partition function and free energy are 

Zo =< exp{-Fo{B)/T) >b, 



"'^private communication to one of the authors (Yu. A.), December 1993. 
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Fo{B)/T = hindetW - lndet{-D\Bj) = (3) 
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where W = —D^{B) — 2gF and D^{B) is the inverse gluon and ghost propagator, respectively, 
is a regularizing factor [12]. 
The ghost propagator can be written as [11,12] 



; =< x\ y^"dse^^'(^)b >= f ds{Dz):^e-''^x,y) (4) 



where standard notations [14,15] are used 

^ ^ aJo ^^^^^ ' ^ Pexpig J B^^dz^, 

and a winding path integral is introduced [12] 

with /3 = 1/T. For the gluon propagator an analogous expression holds true, except that in 
(4) one should insert the gluon spin factor Ppexp2gF into $(a;, y). For a quark propagator the 
sum over windings in (5) acquires a factor (—1)" and the quark spin factor is expga^^F^^i, [12]. 

We are now in a position to make the expansion of Z and F in powers of ga^ (i.e. per- 
turbative expansion in as), and the leading-nonperturbative term Zq,Fq - can be represented 
as a sum of contributions with different behavior of which we systematically will keep the 
leading terms 0(iV2), O(iVc) and O(iVO). 

To describe the temperature of dcconfincmcnt phase transition one should specify the phases 
and compute the associated free energy. For the confining phase to lowest order in Q;s,the free 
energy is given by Eq.(3) plus contribution of the energy density e at zero temperature 

F{1) = eV, - l^V.T^ -TY. A^^^^e--^f- + 0(1^ (6) 

where e is defined by the scale anomaly [18] 



and the next terms in (6) correspond to the contribution of mesons (we keep only pion gas) 
and glueballs. Note that s = 0{N^) while two other terms in (6) are 0{N^). 

For the high temperature or second phase we make an assumption that there all color 
magnetic field correlators are the same as in the first phase, while all color electric fields vanish. 
Since at T = color magnetic correlators (CMC) and color electric correlators (CEC) are equal 
due to the Euclidean 0(4) invariance, one has 

< (f;,(s))2 >=< (f;j2 >e, + < (f;j2 >^„,„; < f' f' (8) 

The string tension a which characterizes confinement is due to the electric fields [20,21], 
e.g. in the plane (14) 

2 

(7 = (7^; = Y / jd''x< Ei{x)^x, 0)Ei{0)^0, x) > +... (9) 
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where dots imply higher order terms in Ei. 

The vanishing of as impUes that gluons and quarks are hberated,and this will contribute 
to the free energy in the dcconfined phase via their closed loop terms (3), (4) with all possible 
windings. The CMC enter via the perimeter contribution < ^{x,x) >= fl (see (4)). As a result 
one has for the second high-temperature phase (cf. [12]) 

^(2) = leVs - (iV,2 - m^Q, ^^N,V,T%n, + 0(iV°) (10) 

where and fig are perimeter terms for quarks and gluons, respectively, the latter was esti- 
mated in [22] from the adjoint Polyakov line. In what follows we replace Q by one for simplicity. 

Comparing (6) and (10), and setting F{1) = F{2) at T = Tc, one finds to order O(A^c)) 
disregarding all meson and glueball contributions 

n /V cxs<F^> \ ^'^ 

T = I 3 327r 



For the standard value oi G2 = ^ < F'^ >= 0.012^6^ [19] (note that for n/ = one should 
use an approximately 3 times larger value of G2 [19]) one has for SU{3) and different values 
of Uf — 0, 2, 4, respectively Tc — 240, 150, 134. This should be compared with lattice data [10] 
Tc{lattice) = 240, 146, 131. The agreement is quite good. Note that at large Nc, one has Tc = 
0{Nj^), i.e. the resulting value of Tc does not depend on Nc in this limit. Hadron contributions 
to Tc are 0(A''~^) and therefore suppressed if Tc is below the Hagedorn temperature, as typically 
happens [23]. 

3. In this section we derive the area law for spatial Wilson loops, expressing spatial string 
tension in terms of CMC. 

To this end we write < W{C) > for any loop as [21] 

2 

< W{C) >= exp[-^- I da^^dapxW) < F^,{u)^{u,u')Fpx{u')^{u\u) > (12) 

-|- higher order cumulants] 

For temporal Wilson loops, in the plane i4,i = 1,2,3, only color electric fields Ei = Fi4 
enter in (12), while for spatial ones in the plane ik;i,k — 1,2,3 there appear color magnetic 
fields Bi — \eikiFki] in the standard way [21], one obtains the area law for large Wilson loops 
of size L, L » T^^^ {T^rn) jg ^j^g magnetic correlation length) 

< W{C) >spatial^ exp{-(JsS^i^) (13) 

where the spatial string tension is [18] 

2 

a^ = ^ j (fx<^ B^{x)^{x, 0)B„(0)$(0, x) > +0(< B^ >) (14) 

and n is the component normal to the plane of the contour, while the last term in (14) denotes 
the contribution of the fourth and higher order cumulants. On general grounds, one can write 
for the integrand in (14) 

< Bi{x)^{x, 0)S,(0)$(0, x) 5i,(L'^(x) + (x) + x^^) - ^i^3^^ (15) 
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and only the term D^{x) enters in (14) [18] 



2 

(Js = ^j (fxD^ix) + 0(< >) (16) 

This holds similarly for the temporal Wilson loop in the plane i4,so that one has the area law 
for T < Tc with temporal string tension 

crE = Yl d^xD'^ix) + 0(< >) (17) 

For T = due to the 0(4) invariance, CEC and CMC coincide and (Te = as- For T > Tc in the 
phase (2) CEC vanish, while CMC change on the scale of the dilaton mass ~ IGey. Therefore 
one expects that ag stays intact till the onset of the dimensional reduction mechanism. This 
expectation is confirmed by the lattice simulation - ag stays constant up to T 1.4Tc [1,2]. 
Recent lattice data [2] show an increase of CTg at T Ri 2Tc for SU(2), which could imply the 
early onset of dimensional reduction. 

4. In this section we consider the qq and gg Green functions G{x, y) at T > and derive 
corresponding Hamiltonians for evolution in the spatial direction. We start with the Feynman- 
Schwinger representation [24] for G{x, y), where for simplicity we neglect spin interaction terms 

/•oo roo 

Gix,y)^J^ dsj^ dse-''-''{Dz)^^{Dz):^<W{C)> (18) 

Here K and {Dz)^y are defined in (5) and < W{C) > in (12), where the contour C is formed 
by paths z{t), z{t) and t = x — y is for definiteness taken along the 3'rd axis. Since by 
definition at T > Tc electric correlators are zero, only elements da^^, in (12) in planes 12,13 and 
23 contribute. As a result one obtains for < W{C) > the form (13) with 

Smin^ dr d-f^Jwlw'^ - {WiW'i)^ (19) 

where only spatial components Wi,i = 1, 2, 3 enter 

Wj(T,7) = 2;j(r)7 + 2;i(r)(l -7),Wi = — , w^ = — (20) 

The form (20) is equivalent to that used before in [14,15] but with = 0. 

As a next step one can introduce the "dynamical masses" fi,fi similarly to [14,15]. We are 
looking for the "cm" Hamiltonian which corresponds to the hyperplane where z^ = ^3. Now 
the role of evolution parameter (time) is played hy Z3 = Z3 = t with < r < t, and we define 
transverse vectors z±{zi, Z2), z±{zi, Z2) and 2:4 (r), Z4{t). 

dzs dr „ dz-i dr ,^ . 



Then K,K m (18) assumes the form 

^^2Jo "^^^W) + ^^^^^^ + + ^^'^^ ^^^^ 
and the same for K with additional bars over n, z±,Z4. 
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Performing the transformation in the functional integral (18) dsDz^{T) — > D^, dsDz^{T) — > 
Dji, one has 

G{x,y) = j DiiDjlDz±Dz4Dz4)^y(Dz4)^yexp(-A) (23) 

with the action 

A^K + K + aS^in (24) 

Note that z^, z^ are not governed by NP dynamics and enter A only kinematically (through 
K,K), and hence can be easily integrated out in (23) using Eq.(5) for the 4-th components - 
with {x — y)4 = (see [18] for more detail). 

One can now proceed as was done in [15,16], i.e. one introduces auxiliary functions z/(r, 7), r;(r, 7), 

— # — * 

defines center-of-mass and relative coordinates R±,f± = r, and finally integrates out R± and 
r]{T,^). The only difference from [15,16] is that now Z4,Zi do not participate in all those 
transformations. As a result one obtains 



G{x,y) ~ / DiyDi^DjlDre-^^'^ ^ g'^^i^^ (25) 

n,n2 



here 

p. .) = i /' drl^i^ + + + p + .V r'-l+r .d,+ (26) 

2 Jo u a Jo 1/ Jo 



Ml - C)2 + + /o' ci7(7 - V 



An.n.--^W ldri^ + -^h (27) 

b{n) = 2n for bosons and 2n + 1 for quarks. We also have introduced the radial momentum 
Pr;the angular momentum L 



IjL + jl 



fr^'^-i-^f'^: L-r-xp (2J 

and 



+ /i + J udj J udj 

Let us define the Hamiltonian H for the given action A = A^^^ + A^^^j (25), integrating 
over Diy, Dfj,, Dft around the extremum of A ( this is an exact procedure in the limit t — > 00). 
For the extremal values of the auxihary fields, one has 

^ + ml + {b{n,)7TTr ^ ^ _ l{l + I) (l-CY _ 1 
/i2(r) al iJ? 

f + ml + {h{n,)7:Tf _ _ + e _ 1 . .... 



z/2(r, 7) 03 



where 03 = /i(l — Q + fi^ + / ^7(7 — C) and C is defined by eq.(29). 

After the substitution of these extremal values into the path integral Hamiltonian 

G{x,y) =< x\ ^ e"-^"i"2*|7/ > (31) 

nin2 
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one has to construct (performing proper Weil ordering [28]) the operator Hamiltonian acting 
on the wave functions. 

Consider for simphcity the case L — O.Then from (30), one obtains 



Hn,n, = ^p' + m\ + (6(ni)7rT)2 + ^p' + ml + (6(n2)7rT)2 + ar (32) 

Here p = and f = f± is a 2d vector, f = (ri,r2); mi,m2 - current masses of quark and 
antiquark and a = as, for the gg system mi = 1712 = and a = Oadj = fcsfor SU(3) [25]. 
Eigenvalues and eigenf unctions of i^mnj 

Hnin2^{r) = M{ni,n2)(p{r) (33) 

define the so-called screening masses and corresponding wave-functions, which have been mea- 
sured in lattice calculations [3-6]. 

The lowest mass sector for mesons is given by i/ooi where for n-i = n2 = one has b{ni) = 
b{n2) = 1 in (32). For light quarks one can put mi = m2 = and expand at large T square 
roots in (32) to obtain 

-2 

//oo~27rT+4^ + (7,r, Mqo ~ 27rr + £(r) (34) 

TTI 

where e{T) — ^^j^^tjtO' 0, ~ 1.74 is the eigenvalue of the 2-D dimensionless Schrodinger 
equation. 

Assuming the parametrization ^ya^ = cg^(T)T with c being numerical constant differ- 
ent for SU{2) and SU{3) [2] and the scahng behavior of g^{T), one has Mqo ^ 2t:T + 
0{{\n{T / Kt))~'^^^T) tending to twice the lowest Matsubara frequency,as was suggested before 
in [26]. (This limit corresponds to the free quarks,propagating perturbativly in the space-time 
with the imposed antipcriodic boundary conditions along the 4th axis). Eq. (34) coincides 
with that proposed in [17], where also numerical study was done of Mqo and 9300 (?")• Our calcu- 
lations of Eq. (33-34) for SU{2>) with c = 0.586 [2] agree with [17] and are presented in Fig.l 
together with lattice calculations of </?(r)oo for p-meson [5]. The values of Moo(T) found on the 
lattice [3,4] are compared with our results in Fig. 2. Note, that our Moo (34) docs not contain 
perimeter corrections which are significant. Therefore one has to add the meson constant to 
Moo to compare with lattice data. We disregard the spin-dependent and one-gluon-exchange 
(OGE) interaction here for lack of space and will treat this subject elsewhere [18]. It is known 
[17], however, that OGE is not very important at around T ^ 2Tc. 

We note that the lowest meson screening mass appears also in Yukawa type exchanges 
between quark lines and can be compared with the corresponding masses in Yukawa potentials. 

5. For the gg system, the lowest mass sector is given by h{n = 0) = 0, and one has from 
(32) 

Hm{gg) = 2]p1 Gadjr (35) 

Note that T does not enter the kinetic terms of (35). 

To calculate with (35) one can use the approximation in (26) of r-independent // [14], which 
leads to the operator (// = p) 

h{n) ^ — + Uadjr (36) 
A* 

The eigenvalue E^jj) of h{ii) should be minimized with respect to ji and the result E{fio) is 
known to yield the eigenvalue of (35) within a few percent accuracy [27]. The values E{iiq) ^ 
Mgg thus found are presented in Fig. 2. The corresponding wave functions (/^oo(r) are given in 
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Fig.l. These data can be compared with the glueball screening masses, found on the lattice in 
[6]. 

Another point of comparison is the Polyakov hne correlator 
with 

ViR) = , niR) = -^trPexpitg f A,dz. 

K" I\c JO 

It is easy to understand that M = Mgg and one can compare our results for the GSM Mgg(T) 
with the corresponding lattice data for M in Fig. 2. 

Taking again into account an unknown constant perimeter correction to our values of Mgg, 
one can see a reasonable qualitative behavior. We refer the reader to another publication for 
more detail [18]. 

In addition to the purely nonperturbative source of GSM described above there is a com- 
peting mechanism - the perturbative formation of the electric Debye screening mass mei{T) = 
gT{^ + ^) for each gluon with g{T) given by the temperature dependence of ^JajT) = 
cg'^{T)T [2]. Therefore for large T, where rriei is essential, one should rather use instead of (35) 
another Hamiltonian, which is obtained from (32) replacing mi = m2 = rrieiiT) 



Hooigg) = 2^f + ml{T) + <j„,,(T)r (38) 

It turns out that up to temperatures T < 2Tc, one can consider the effect of rrieiiT) as a small 
correction to the eigenvalue of (35) E 

Mgg^E + 5E^ 4(^)3/^ . ^c/(T)r + (39) 

where a ^ 1.74 and this justifies our nonperturbative treatment of the GSM in the temperature 
region concerned. 

We note that the transition between these two regimes (with the dominance of nonpertur- 
bative and then perturbative dynamics) is smooth, with the GSM tending to twice the Debye 
mass at large T [18]. 

6. In conclusion, we have pursued a selfconsistent analysis of nonperturbative dynamics 
at T > Tc, including the calculation of Tc, spatial area law and meson and glueball screening 
masses. We stress that in both meson and glueball cases, the nonperturbative mechanism 
plays an essential role in the formation of the screening masses in the temperature region 
< T < 2-=-3Tc, while at large T, the lowest screening masses tend to the corresponding limits, 
governed by perturbative dynamics. We obtain in all cases a reasonable agreement with known 
lattice data supporting our picture of magnetic confinement above Tc, while there is electric 
confinement below Tc. An additional piece of data is necessary to confirm our picture, e.g. 
lattice calculation of the screening glueball wave function to compare with our prediction. 

One of us (E.G.) is grateful to B.L. loffe and V. Eletskii for useful discussions. The authors 
are grateful to D. Melikhov for help in computations. This work was supported by the Russian 
Fund for Fundamental Research, grant 93-02-14937. 
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Figure 1: Wave functions for the p-meson (solid line) and glueball (dashed line) vs r/a (a=0.23 
fm) for T = 210 MeV. The data from ref.[5] correspond to T ^ 210 MeV. 




Figure 2: The screening masses for mesons (solid line) and glueballs (dashed line) as a function 
of temperature. The data from refs. [3,4] (squares) and [6] (triangle). 
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